regioneReloaded: evaluating the association of multiple genomic region sets.

MOTIVATION: Next-generation sequencing methods continue improving the annotation of genomes in part by determining the distribution of features such as epigenetic marks. Evaluating and interpreting the association between genomic regions and their features has become a common and challenging analysis in genomic and epigenomic studies. RESULTS: With regioneR we provided an R package allowing to assess the statistical significance of pairwise associations between genomic region sets using permutation tests. We now present the R package regioneReloaded that builds upon regioneR's statistical foundation and extends the functionality for the simultaneous analysis and visualization of the associations between multiple genomic region sets. Thus, we provide a novel discovery tool for the identification of significant associations that warrant to be tested for functional interdependence. AVAILABILITY AND IMPLEMENTATION: regioneReloaded is an R package released under an Artistic-2.0 License. The source code and documentation are freely available through Bioconductor: http://www.bioconductor.org/packages/regioneReloaded.

The consequences of cohesin mutations in myeloid malignancies.

Recurrent somatic mutations in the genes encoding the chromatin-regulatory cohesin complex and its modulators occur in a wide range of human malignancies including a high frequency in myeloid neoplasms. The cohesin complex has a ring-like structure which can enclose two strands of DNA. A first function for the complex was described in sister chromatid cohesion during metaphase avoiding defects in chromosome segregation. Later studies identified additional functions of the cohesin complex functions in DNA replication, DNA damage response, 3D genome organisation, and transcriptional regulation through chromatin looping. In this review, we will focus on STAG2 which is the most frequently mutated cohesin subunit in myeloid malignancies. STAG2 loss of function mutations are not associated with chromosomal aneuploidies or genomic instability. We hypothesize that this points to changes in gene expression as disease-promoting mechanism and summarize the current state of knowledge on affected genes and pathways. Finally, we discuss potential strategies for targeting cohesion-deficient disease cells.

The histone chaperone ANP32B regulates chromatin incorporation of the atypical human histone variant macroH2A.

All vertebrate genomes encode for three large histone H2A variants that have an additional metabolite-binding globular macrodomain module, macroH2A. MacroH2A variants impact heterochromatin organization and transcription regulation and establish a barrier for cellular reprogramming. However, the mechanisms of how macroH2A is incorporated into chromatin and the identity of any chaperones required for histone deposition remain elusive. Here, we develop a split-GFP-based assay for chromatin incorporation and use it to conduct a genome-wide mutagenesis screen in haploid human cells to identify proteins that regulate macroH2A dynamics. We show that the histone chaperone ANP32B is a regulator of macroH2A deposition. ANP32B associates with macroH2A in cells and in vitro binds to histones with low nanomolar affinity. In vitro nucleosome assembly assays show that ANP32B stimulates deposition of macroH2A-H2B and not of H2A-H2B onto tetrasomes. In cells, depletion of ANP32B strongly affects global macroH2A chromatin incorporation, revealing ANP32B as a macroH2A histone chaperone.

Histone Variants and Their Chaperones in Hematological Malignancies.

Epigenetic regulation occurs on the level of compacting DNA into chromatin. The functional unit of chromatin is the nucleosome, which consists of DNA wrapped around a core of histone proteins. While canonical histone proteins are incorporated into chromatin through a replication-coupled process, structural variants of histones, commonly named histone variants, are deposited into chromatin in a replication-independent manner. Specific chaperones and chromatin remodelers mediate the locus-specific deposition of histone variants. Although histone variants comprise one of the least understood layers of epigenetic regulation, it has been proposed that they play an essential role in directly regulating gene expression in health and disease. Here, we review the emerging evidence suggesting that histone variants have a role at different stages of hematopoiesis, with a particular focus on the histone variants H2A, H3, and H1. Moreover, we discuss the current knowledge on how the dysregulation of histone variants can contribute to hematopoietic malignancies.

Evolution, structure and function of divergent macroH2A1 splice isoforms.

The replacement of replication-coupled histones with non-canonical histone variants provides chromatin with additional properties and contributes to the plasticity of the epigenome. MacroH2A histone variants are counterparts of the replication-coupled histone H2A. They are characterized by a unique tripartite structure, consisting of a histone fold, an unstructured linker, and a globular macrodomain. MacroH2A1.1 and macroH2A1.2 are the result of alternative splicing of the MACROH2A1 gene and can have opposing biological functions. Here, we discuss the structural differences between the macrodomains of the two isoforms, resulting in differential ligand binding. We further discuss how this modulates gene regulation by the two isoforms, in cases resulting in opposing role of macroH2A1.1 and macroH2A1.2 in development and differentiation. Finally, we share recent insight in the evolution of macroH2As. Taken together, in this review, we aim to discuss in unprecedented detail distinct properties and functions of the fascinating macroH2A1 splice isoforms.

Different Gene Sets Are Associated With Azacitidine Response In Vitro Versus in Myelodysplastic Syndrome Patients.

Myelodysplastic syndromes (MDS) are a heterogeneous group of hematopoietic disorders characterized by dysplasia, ineffective hematopoiesis, and predisposition to secondary acute myeloid leukemias (sAML). Azacitidine (AZA) is the standard care for high-risk MDS patients not eligible for allogenic bone marrow transplantation. However, only half of the patients respond to AZA and eventually all patients relapse. Response-predicting biomarkers and combinatorial drugs targets enhancing therapy response and its duration are needed. Here, we have taken a dual approach. First, we have evaluated genes encoding chromatin regulators for their capacity to modulate AZA response. We were able to validate several genes, whose genetic inhibition affected the cellular AZA response, including 4 genes encoding components of Imitation SWItch chromatin remodeling complex pointing toward a specific function and co-vulnerability. Second, we have used a classical cohort analysis approach measuring the expression of a gene panel in bone marrow samples from 36 MDS patients subsequently receiving AZA. The gene panel included the identified AZA modulators, genes known to be involved in AZA metabolism and previously identified candidate modulators. In addition to confirming a number of previously made observations, we were able to identify several new associations, such as NSUN3 that correlated with increased overall survival. Taken together, we have identified a number of genes associated with AZA response in vitro and in patients. These groups of genes are largely nonoverlapping suggesting that different gene sets need to be exploited for the development of combinatorial drug targets and response-predicting biomarkers.

The high mobility group protein HMG20A cooperates with the histone reader PHF14 to modulate TGFß and Hippo pathways.

High mobility group (HMG) proteins are chromatin regulators with essential functions in development, cell differentiation and cell proliferation. The protein HMG20A is predicted by the AlphaFold2 software to contain three distinct structural elements, which we have functionally characterized: i) an amino-terminal, intrinsically disordered domain with transactivation activity; ii) an HMG box with higher binding affinity for double-stranded, four-way-junction DNA than for linear DNA; and iii) a long coiled-coil domain. Our proteomic study followed by a deletion analysis and structural modeling demonstrates that HMG20A forms a complex with the histone reader PHF14, via the establishment of a two-stranded alpha-helical coiled-coil structure. siRNA-mediated knockdown of either PHF14 or HMG20A in MDA-MB-231 cells causes similar defects in cell migration, invasion and homotypic cell-cell adhesion ability, but neither affects proliferation. Transcriptomic analyses demonstrate that PHF14 and HMG20A share a large subset of targets. We show that the PHF14-HMG20A complex modulates the Hippo pathway through a direct interaction with the TEAD1 transcription factor. PHF14 or HMG20A deficiency increases epithelial markers, including E-cadherin and the epithelial master regulator TP63 and impaired normal TGFß-trigged epithelial-to-mesenchymal transition. Taken together, these data indicate that PHF14 and HMG20A cooperate in regulating several pathways involved in epithelial-mesenchymal plasticity.

MacroH2As regulate enhancer-promoter contacts affecting enhancer activity and sensitivity to inflammatory cytokines.

MacroH2A histone variants have a function in gene regulation that is poorly understood at the molecular level. We report that macroH2A1.2 and macroH2A2 modulate the transcriptional ground state of cancer cells and how they respond to inflammatory cytokines. Removal of macroH2A1.2 and macroH2A2 in hepatoblastoma cells affects the contact frequency of promoters and distal enhancers coinciding with changes in enhancer activity or preceding them in response to the cytokine tumor necrosis factor alpha. Although macroH2As regulate genes in both directions, they globally facilitate the nuclear factor κB (NF-κB)-mediated response. In contrast, macroH2As suppress the response to the pro-inflammatory cytokine interferon gamma. MacroH2A2 has a stronger contribution to gene repression than macroH2A1.2. Taken together, our results suggest that macroH2As have a role in regulating the response of cancer cells to inflammatory signals on the level of chromatin structure. This is likely relevant for the interaction of cancer cells with immune cells of their microenvironment.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity.

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

PLCG1 is required for AML1-ETO leukemia stem cell self-renewal.

In an effort to identify novel drugs targeting fusion-oncogene-induced acute myeloid leukemia (AML), we performed high-resolution proteomic analysis. In AML1-ETO (AE)-driven AML, we uncovered a deregulation of phospholipase C (PLC) signaling. We identified PLCgamma 1 (PLCG1) as a specific target of the AE fusion protein that is induced after AE binding to intergenic regulatory DNA elements. Genetic inactivation of PLCG1 in murine and human AML inhibited AML1-ETO dependent self-renewal programs, leukemic proliferation, and leukemia maintenance in vivo. In contrast, PLCG1 was dispensable for normal hematopoietic stem and progenitor cell function. These findings are extended to and confirmed by pharmacologic perturbation of Ca++-signaling in AML1-ETO AML cells, indicating that the PLCG1 pathway poses an important therapeutic target for AML1-ETO+ leukemic stem cells.

Evolution of a histone variant involved in compartmental regulation of NAD metabolism.

NAD metabolism is essential for all forms of life. Compartmental regulation of NAD+ consumption, especially between the nucleus and the mitochondria, is required for energy homeostasis. However, how compartmental regulation evolved remains unclear. In the present study, we investigated the evolution of the macrodomain-containing histone variant macroH2A1.1, an integral chromatin component that limits nuclear NAD+ consumption by inhibiting poly(ADP-ribose) polymerase 1 in vertebrate cells. We found that macroH2A originated in premetazoan protists. The crystal structure of the macroH2A macrodomain from the protist Capsaspora owczarzaki allowed us to identify highly conserved principles of ligand binding and pinpoint key residue substitutions, selected for during the evolution of the vertebrate stem lineage. Metabolic characterization of the Capsaspora lifecycle suggested that the metabolic function of macroH2A was associated with nonproliferative stages. Taken together, we provide insight into the evolution of a chromatin element involved in compartmental NAD regulation, relevant for understanding its metabolism and potential therapeutic applications.

Divergent leukaemia subclones as cellular models for testing vulnerabilities associated with gains in chromosomes 7, 8 or 18.

Haematopoietic malignancies are frequently characterized by karyotypic abnormalities. The development of targeted drugs has been pioneered with compounds against gene products of fusion genes caused by chromosomal translocations. While polysomies are equally frequent as translocations, for many of them we are lacking therapeutic approaches aimed at synthetic lethality. Here, we report two new cell lines, named MBU-7 and MBU-8, that differ in complete trisomy of chromosome18, a partial trisomy of chromosome 7 and a tetrasomy of the p-arm of chromosome 8, but otherwise share the same mutational pattern and complex karyotype. Both cell lines are divergent clones of U-937 cells and have the morphology and immunoprofile of monocytic cells. The distinct karyotypic differences between MBU-7 and MBU-8 are associated with a difference in the specific response to nucleoside analogues. Taken together, we propose the MBU-7 and MBU-8 cell lines described here as suitable in vitro models for screening and testing vulnerabilities that are associated with the disease-relevant polysomies of chromosome 7, 8 and 18.

Inhibition of CBP synergizes with the RNA-dependent mechanisms of Azacitidine by limiting protein synthesis.

The nucleotide analogue azacitidine (AZA) is currently the best treatment option for patients with high-risk myelodysplastic syndromes (MDS). However, only half of treated patients respond and of these almost all eventually relapse. New treatment options are urgently needed to improve the clinical management of these patients. Here, we perform a loss-of-function shRNA screen and identify the histone acetyl transferase and transcriptional co-activator, CREB binding protein (CBP), as a major regulator of AZA sensitivity. Compounds inhibiting the activity of CBP and the closely related p300 synergistically reduce viability of MDS-derived AML cell lines when combined with AZA. Importantly, this effect is specific for the RNA-dependent functions of AZA and not observed with the related compound decitabine that is only incorporated into DNA. The identification of immediate target genes leads us to the unexpected finding that the effect of CBP/p300 inhibition is mediated by globally down regulating protein synthesis.

The Histone Variant MacroH2A1 Impacts Circadian Gene Expression and Cell Phenotype in an In Vitro Model of Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death worldwide. A foremost risk factor for HCC is obesity/metabolic syndrome-related non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH), which is prompted by remarkable changes in transcription patterns of genes enriching metabolic, immune/inflammatory, and circadian pathways. Epigenetic mechanisms play a role in NAFLD-associated HCC, and macroH2A1, a variant of histone H2A, is involved in the pathogenesis modulating the expression of oncogenes and/or tumor suppressor genes and interacting with SIRT1, which crucially impacts the circadian clock circuitry. Hence, we aimed to appraise if and how macroH2A1 regulated the expression patterns of circadian genes in the setting of NAFLD-associated HCC. We took advantage of an in vitro model of liver cancer represented by HepG2 (human hepatocarcinoma) cells stably knocked down for macroH2A1 and conducted whole transcriptome profiling and deep phenotyping analysis. We found up-regulation of PER1 along with several deregulated circadian genes, enriching several important pathways and functions related to cancer onset and progression, such as epithelial-to-mesenchymal transition, cell cycle deregulation, and DNA damage. PER1 silencing partially mitigated the malignant phenotype induced by the loss of macroH2A1 in HCC cells. In conclusion, our findings suggest a modulatory role for the core circadian protein PER1 in liver carcinogenesis in the context of a lack of the macroH2A1 epigenetic and transcriptional landscape.

The Role of MacroH2A Histone Variants in Cancer.

The epigenome regulates gene expression and provides a molecular memory of cellular events. A growing body of evidence has highlighted the importance of epigenetic regulation in physiological tissue homeostasis and malignant transformation. Among epigenetic mechanisms, the replacement of replication-coupled histones with histone variants is the least understood. Due to differences in protein sequence and genomic distribution, histone variants contribute to the plasticity of the epigenome. Here, we focus on the family of macroH2A histone variants that are particular in having a tripartite structure consisting of a histone fold, an intrinsically disordered linker and a globular macrodomain. We discuss how these domains mediate different molecular functions related to chromatin architecture, transcription and DNA repair. Dysregulated expression of macroH2A histone variants has been observed in different subtypes of cancer and has variable prognostic impact, depending on cellular context and molecular background. We aim to provide a concise review regarding the context- and isoform-dependent contributions of macroH2A histone variants to cancer development and progression.

Poly(ADP-ribose) binding and macroH2A mediate recruitment and functions of KDM5A at DNA lesions.

The histone demethylase KDM5A erases histone H3 lysine 4 methylation, which is involved in transcription and DNA damage responses (DDRs). While DDR functions of KDM5A have been identified, how KDM5A recognizes DNA lesion sites within chromatin is unknown. Here, we identify two factors that act upstream of KDM5A to promote its association with DNA damage sites. We have identified a noncanonical poly(ADP-ribose) (PAR)-binding region unique to KDM5A. Loss of the PAR-binding region or treatment with PAR polymerase (PARP) inhibitors (PARPi's) blocks KDM5A-PAR interactions and DNA repair functions of KDM5A. The histone variant macroH2A1.2 is also specifically required for KDM5A recruitment and function at DNA damage sites, including homology-directed repair of DNA double-strand breaks and repression of transcription at DNA breaks. Overall, this work reveals the importance of PAR binding and macroH2A1.2 in KDM5A recognition of DNA lesion sites that drive transcriptional and repair activities at DNA breaks within chromatin that are essential for maintaining genome integrity.

Disruption of paternal circadian rhythm affects metabolic health in male offspring via nongerm cell factors.

Circadian rhythm synchronizes each body function with the environment and regulates physiology. Disruption of normal circadian rhythm alters organismal physiology and increases disease risk. Recent epidemiological data and studies in model organisms have shown that maternal circadian disruption is important for offspring health and adult phenotypes. Less is known about the role of paternal circadian rhythm for offspring health. Here, we disrupted circadian rhythm in male mice by night-restricted feeding and showed that paternal circadian disruption at conception is important for offspring feeding behavior, metabolic health, and oscillatory transcription. Mechanistically, our data suggest that the effect of paternal circadian disruption is not transferred to the offspring via the germ cells but initiated by corticosterone-based parental communication at conception and programmed during in utero development through a state of fetal growth restriction. These findings indicate paternal circadian health at conception as a newly identified determinant of offspring phenotypes.

Epigenetics in a Spectrum of Myeloid Diseases and Its Exploitation for Therapy.

Mutations in genes encoding chromatin regulators are early events contributing to developing asymptomatic clonal hematopoiesis of indeterminate potential and its frequent progression to myeloid diseases with increasing severity. We focus on the subset of myeloid diseases encompassing myelodysplastic syndromes and their transformation to secondary acute myeloid leukemia. We introduce the major concepts of chromatin regulation that provide the basis of epigenetic regulation. In greater detail, we discuss those chromatin regulators that are frequently mutated in myelodysplastic syndromes. We discuss their role in the epigenetic regulation of normal hematopoiesis and the consequence of their mutation. Finally, we provide an update on the drugs interfering with chromatin regulation approved or in development for myelodysplastic syndromes and acute myeloid leukemia.

Histone Modifications and Their Targeting in Lymphoid Malignancies.

In a wide range of lymphoid neoplasms, the process of malignant transformation is associated with somatic mutations in B cells that affect the epigenetic machinery. Consequential alterations in histone modifications contribute to disease-specific changes in the transcriptional program. Affected genes commonly play important roles in cell cycle regulation, apoptosis-inducing signal transduction, and DNA damage response, thus facilitating the emergence of malignant traits that impair immune surveillance and favor the emergence of different B-cell lymphoma subtypes. In the last two decades, the field has made a major effort to develop therapies that target these epigenetic alterations. In this review, we discuss which epigenetic alterations occur in B-cell non-Hodgkin lymphoma. Furthermore, we aim to present in a close to comprehensive manner the current state-of-the-art in the preclinical and clinical development of epigenetic drugs. We focus on therapeutic strategies interfering with histone methylation and acetylation as these are most advanced in being deployed from the bench-to-bedside and have the greatest potential to improve the prognosis of lymphoma patients.